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ABSTRACT 

To determine the equation of state of the universe, we propose to use a new independent 
variable R = {Ho/c){dL{z)/{l + z)), where Hq and dL{z) are the present Hubble parameter 
and the luminosity distance, respectively. For the flat universe suggested from the observation 
of the anisotropy of cosmic microwave background, the density and the pressure are expressed 

as p/po — 4:{df / dR)"^ / and p/po = —A/'d{d^f/dR^)/f^ where po is the present density and 
f{R) = + z{R). In (R, f) plane the sign as well as the strength of the pressure is in 

proportion to the curvature of the curve f{R)- We propose to adopt a Pade-like expression of 

N 

f{R) = ^/y/u with u = 1 + ^ For flat A model the expansion up to = 7 has at 

n=l 

most an error < 0.2% for z < 1.7 and any value of A. We also propose a general method to 
determine the equation of state of the universe which has N — 1 free parameters. If the number 
of parameters are smaller than N — 1, there is a consistency check of the equation of state so 
that we may confirm or refute each model. 

Subject headings: cosmology: theory, distance scale, dark matter 
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Introduction 



Recent measurements of the luminosity distance dL{z) using Type la supernovae (Schmidt et al 



1998, Riess et al. 1998, Pcrlmutter et al. 1999) suggest that accurate d^lz) may be obtained in the near 



future. Especially SNAP (Supernova/ Acceleration Probe) will give us the luminosity distance of ~2000 
Type la supernovae with an accuracy of a few % up to z ^ 1.7 every year. On the other hand from the 
observation of the first Doppler peak of the anisotropy of CMB, it is now suggested that the universe is 
flat ( |de Bernardis et al. 200C , Lange et al. 200C ), which may be proved in future by MAP and Planck. If 
the flat universe is the case, the density p{z) and the pressure p{z) can be determined only from dL{z) in 



principle (Nakamura & Chiba 1999) so that the equation of the state of the universe is uniquely determined. 
If not, the determination of the present curvature of the universe and the determination of p{z) and p{z) 
will be coupled in general ( Nakamura fc Chiba 1999| ). 



Now let us assume that the universe is flat. Even in this case at least two problems exist: 1) How 
to express the continuous function dL{z) which is accurate enough from z = to z ^ 1.7 using several 
free parameters; 2) From dL{z) how to obtain accurate p{z) and p{z), that is, the equation of state of the 



universe QStarobinsky 1998| , [Huterer fc Turner 1999| , [Nakamura fc Chiba 1999| , [Saini et al. 2000| , phiba fc 
Nakamura 2"000| ). 



In this paper, we propose to use a new independent variable R = (Hq / c){dL{z) / (1 + z)) instead 
of z, where Hq is the present Hubble parameter. We show that p{R)/po = 4:{df / dR)'^ / and 
p{R)/Pq — ~4:/3{d^f/dR'^)/f^ where po is the present density and f{R) — l/yjl + z{R). This means that 
the pressure is in proportion to the curvature of the curve f{R). For an accurate expression of f{R) we 

N 

propose a Pade-like form of f{R) = l/\Ai with m = 1 + ^ UnR^. For flat A model, the expansion up to 

■11=1 

N — 7 has at most an error < 0.2% for z < 1.7 and all value of A. We also propose a general method to 
determine the equation of state of the universe which has < — 1 free parameters. 



2. New Variables 

The luminosity distance dL{z) is given by dL{z) ~ ao(l + z)f{x), with x — l/^o dz' / H{z'), where 
H{z) and oq are the Hubble parameter at z and the present scale factor, and /(x) = Xi sinh(x) and sin(x) 
for flat, open and closed universe, respectively. Let us define R = {Ha/ c){di^{z) / {1 + z)), where Hq is the 
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present Hubble parameter. Then p{z) and p{z) are expressed as 

^ - -4r- + (1 + f t1 - f(i + - t + (1 + -):^ f t1 1 ^o^^^o, (1) 



where po is the present density and ilfeo = k/{aQHQ) ( Nakamura fc Chiba 199£ ). Since the fiat universe 
is suggested both from observations (|de Bernardis et al. 200d| , Lange et al. 2000|) and theory (inflation 
paradigm), we consider only — case in this paper. 

Now we adopt R as an independent variable instead of z. Then p{R) and p{R) are expressed 



as p{R)/po = 4:{df/dR)^f-^ andp{R)/po = -4:/3{d^f/dR'^)f-^, where f{R) = 1/^1 + z{R). These 
expressions of p{R) and p{R) have quite interesting physical meanings. The density is in proportion to the 
square of the first derivative of / with respect to the new independent variable R, while the pressure is 
in proportion to the second derivative of /, that is, the curvature. Therefore if the pressure is zero, f{R) 
is the straight line while the negative pressure corresponds to the positive curvature of the curve f{R) in 
{R, f) plane. This is completely in contrast to (z,dL{z)) plane where as far as p > 3p, the curve dL{z) 
has the positive curvature. Therefore in plane it is difficult to distinguish by eye if the pressure 

is negative or not. However in (i?, /) plane it is quite easy to distinguish the sign of the curvature of the 
curve so that the negative pressure can be identified by eye. To demonstrate this we show in Fig.l. the 
data from Supernova Cosmology Project ( Perlmutter et al. 1999 ) in {R, f) plane.0 Although the error bar 



is considerably large, a glance may show that the curvature is positive. 

To analyze the data quantitatively, the first way is to expand / in a power series of 

N 

R as f{R) = 1 + X] fnR", where /„ are constants. We did this expansion with = 4 



for the data from Supernova Cosmology Project ( Perlmutter et al. 1999|) and obtained Q 



/i = -0.5, /2 ^ 0.292to;oi5, /s ^ -0.257 ± 0.034, and /4 = 0.046^°;^?^. Therefore the 
apparent positive curvature by eye is confirmed quantitatively since /2 is positive. ^ 



^ Similar analysis may be possible using data from High-Z Supernova Search (Schmidt et al. 1998). As 
an example we use only the data from Supernova Cosmology Project in this letter. This does not mean that 
other data are not important. 

^In the likelihood analysis, we used the inverse function R = R{z) by solving the quartic equation. 

■^Note here that for the flat A model with — 0.3 and z — 1 which is roughly corresponding to the data 
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To know the convergence property of the expansion of / in a power series of R, let us consider the flat 
A model. In this model, R is expressed as 



i? = 2 / df/^QM + (1 - nM)f^ (2) 



/ 

where ^Im is the present density parameter of the dark matter. For N — 7 the expansion is expressed as 

/ = i-lR + l{i-nM)R'-^{i-nM)R' 

128 zoo 
+ - nM)mM + 168(1 - ^m) + 12r!Af (1 - ^m) + 63(1 - nM)^]R'' 

-T7^(l - ^M)mM + 1176(1 - Qm) + 132r!M(l - ^m) + 1827(1 - (3) 
14336 

However in this expansion the sign of /« changes alternately so that the convergence is extremely slow. The 
accuracy of the expansion is only 16% for Hm = 0.3 and / = 0.6 corresponding to z — 1.77 and R ~ 1.13. 
While for low z, the accuracy is better. For example for Qm = 0.3 and / = 0.7(i? = 0.79, z = 1.04) the 
accuracy is 0.95%. 

The above expansion might not be accurate enough to analyze the data from such as SNAP so that 
we need a Fade-like approximation in a power series expansion of /. We adopt the following Fade- like 
approximation as 

N 

/(i?)-l/V^, U=l + J2nnR". (4) 

n=l 

For the flat A model, R is related to u as 



R = J du/^QMU^ + (1 - ^m). (5) 

Up to A'^ = 7, w is expanded as 

u = i + R + ^nMR^ + lnMR'' + ^nM{2 + 3nM)R^ 

A Z lb 

+^nljR' + ^f]2^(4 + 3nM)R' + YY^f^M(i + Q^ai)R'- (6) 

All u„ are positive definite so that the convergence is very rapid. In reality for — 0.3 and 

/ = 0.6(i? = 1.13, u = 2.77) the relative error of / is 0.059%. In this expansion the error is the largest 



of Supernova Cosmology Project ( Perlmutter et al. 1999D , the expansion of / with A^ = 4 has an accuracy 



of 3.6%. Here the accuracy means the relative error of / computed from the expansion (3) to the given / in 
Eq. (2). 
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for flM = 1- However even in this case, the error is 0.2% for / = 0.6(i? = 0.8, u — 2.77). Therefore 
if accurate vahies of R for various z are obtained observationaUy, we may determine seven parameters 
Un{n— 1, 2, ... 7), which may be enough to express f{R) in < 0.1% accuracy. 



3. Equation of State 

Let us assume that accurate u„ for n = 1, . . . , N are obtained from the analysis of observational data. 
In this section we discuss how to determine and confirm the equation of state from m„. Firstly p{R) and 
p{R) are expressed as p{R)/po — (du/dR)^ and p{R)/po = —(du/dR)'^ + 2u/3{d'^u/dR'^). Secondly the 
square of the sound velocity dp/ dp is expressed as 



dp_ 2 uf " 



(7) 



' dR dW 



3.1. One parameter Equation of State 

The flat A model has only one parameter, r2j\/. This parameter is equal to either 4u2/3 or 2it3 so 
that the identity of Vlj^i = 4u2/3 = 2u3 will be the consistency check of A model (see also [Chiba fc" 



Nakamura 1998). From the data of Supernova Cosmology Project (Perlmutter et al. 199£) we obtained 



^ l^U2 ^ 0.084t°;^^^,U3 = 0.343+°;};^, and = 0.360l^j^ Therefore 4^2/3 = 0.112+°;^°;^ while 
2u3 = 0.686lQ'2g. Note here that for the flat A model with Qm — 0.3 and z = 1 which is roughly 



corresponding to the data of Supernova Cosmology Project (Perlmutter et al. 1999), the expansion of u 
with iV = 4 has an accuracy of 0.19%. There is another null test of the A model for all R. From Eq.(5) we 
can derive 



du/dR = \/Qmu^ + (1 - ^m), 
d^u/dR^ = ^^MU^, 



d^u/dR^ = SnMU^/nMU^ + (1 - ^m)- (8) 

Using these expressions of derivatives, we can easily prove dp/ dp = for all u. Therefore null test of dp/ dp 
in Eq.(7) observationaUy can confirm the A model. At present from the above values of U2 and U3, we have 
dp/ dp at M = 1 as dp/ dp = 3.41^g'2, which means the present data is not accurate enough to confirm the A 
model. 
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3.2. Two pareimeter Equation of State 

The flat w-cosmology is an example of this class. In w-cosmology the universe contains x-matter 



with px = wpx where w is a constant. R is expressed as i? = du/ ^Q.mu^ + (1 ~ ^m)u^^^^'^'> . Then 
derivatives are given by 

du 



^'^MU^ + (1 - i^M)M3(l+'"), 

_ = [Z^j^u + -(1 + w){-iw + 2)(1 - f7M)w'™+ V0mu3 + (1 - 17m)m3(i+-). (9) 



U2 and U3 are given by 



W2 = ^ + 1 , «3 = ^ + 1 • (10) 

Om and w are determined from U2 and U3 as 

(4^2 - 3)2 4^3-2 5 

S2m = 1 7^ 7, = • (11) 

12u3 - 6 - 5(4u2 - 3) ' 4u2-3 3 ^ ' 

Then dp/dp (= w;(l + w){l - ^mW^^"" /{nMU^ + (1 + w){l - ^mW^^"")) should agree with dp/dp in Eq. 
(7) so that we can confirm or refute w-cosmology. 



3.3. General Equation of State 



Let us now consider the general case in which w is a function of u in w-cosmology. R is expressed as 

r du . X r 3(i+w(u)) , 

R= / , g(u) = / ^ ^ " du. (12) 

Ji V^MW^ + (1 - f^M) exp (g'(u)) Ji u 

Then derivatives are given by 

3^ = \/f2Mw3 + (1 - r2M)exp(q(u)), 
oix 

d^u 3_ 9 3,, ^ ,(l+ti;(u)) , , 

= -Qmu^ + o 1 - ^m)- ^ exp {q{u)), 

dR'' 2 2 u 

^ = \Z^MU + -(} ^^^(3ti;(«) + 2) + — )(1 - 17m) exp {q{u))\^J^MV? + (1 - ^m) exp (g(«)). 

U2 and U3 are given by 

3nM , 3(1 + wo)(1-17m) Om , (l-nM)[(l+wo)(3wo + 2)+wi] 
W2 = ^ + ^ , u^ = —^ ^ , (13) 



- 8 - 



where wq = w(l) and wi = dw/du{u — 1). Now we have only two equations for three unknown constants 
^M, W{) and wi. To resolve this one may use the expression for M4. However in the expression of U4, a 
new unknown constant W2 = (Pw / dv?{u ~ 1) appears so that we have to make the closure. One way is 
to determine Q.m from another data such as fii(z) for z > 3 by NGST(Next Generation Space Telescope 



( Efstathiou 1999 )). In this case wq and wi . . . are determined as 



^o-^TT — FTT-l, u'l =4u3-2f7M- (l + ii'o)(3wo + 2), = . . . (14) 

6{l — \Im) 

The other is to assume W2 — 0. Then VLm^wq and wi are determined as a solution to simultaneous non- linear 
equations as 

U2 = li2(f^Af,t«0, Wl), W3 = U3(riAf, Wo,Wi), ^ Ui{Q.M ,Wo,Wi). (15) 

In both cases dp /dp is given by 



dp _ w{u){l + w{u)) + |wi(l - ^m) exp (g(u)) 



(16) 



dp nMU^ + (1 - r!M)i±^ exp {q{u)) ' 

where w(u) = Wq + Wi(u — 1). As before this dp /dp should agree with dp /dp in Eq. (7), which is used 
to confirm or refute each model. In general, if u„ are determined up to iV, the equation of state with 

N-l 

w = Wq + ^ Wn{u — 1)" Can bc determined in principle. 



4. Discussion 

Now if x-matter consists of the scalar field (j) with the potential F (</>), they are related to px and px as 

^(0) = ^{px-px)- 



Using p{R) and p(i?), we have 



If 2 d^^t 

-?,nMU+-—.dR = g{R), (17) 



V8^7o V Mdi?2 



yi<t>) - ^ f 2(^)2 - ^.,u^ - ) - Mi?), (18) 



Vm^ V c?i? 3 di?2 

where 0o is the present value of the scalar field . From Eq. (17) we have R — g^^{4> — <Po)- Then the 
potential is expressed as V{<j)) — h{g^^{(l> — <j)o))- 

In section 3 we assumed that accurate ii„ for n = I, . . . , N are obtained. We here show an example 
of the determination of u„ for n = 1, . . . , 7. We adopt the redshifts of 38 data with z > 0.17 from 
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(Pcrlmuttcr et al. 1999). We also adopt the relative error of R for each data as Xi. Now let us assume 
that our universe obeys the A model with VIm = 0.3. Then we know the theoretical value of i?* for 
each Zi. To simulate the real observation, we set Ri — -R*(l ± SXi) where S is a scale factor. Let us 
assume that an accurate observation gives us 38 luminosity distances with ^ 0.1% accuracy so that 
the scale factor S is chosen to make the relative error of R for 38 data be ~ 0.1%.^ We performed 
the likelihood analysis for this simulated data and obtained U2 = 0.2240, its = 0.1505, U4 = 0.0585, 
W5 — 0.0187,1*6 — 0.0095 and U7 — 0.0053 with = 36.16 for 38-7 d.o.f. while theoretical values are 
U2 = 0.225, U3 = 0.15,1*4 0.054375, 0.016875,1*6 = 0.0068906 and uj = 0.00225. ^From this 

— 4/3ii2 = 0.29866 or flM — 2it3 — 0.3010 is obtained while dp/dp — 0.0052 at it = 1. Assuming the 
more general equation of state with W2 ~ W3 . . . = wi = 0, we have JIm ~ 0.31095,1^0 = —1.01783 and 
wi = —0.04768. This suggests that we may confirm the A model if such accurate luminosity distances are 
available. We show in Fig. 2 the simulated data and the results of the likelihood analysis. Note that error 
bars are extended by a factor 100. The theoretical curve (the dashed line) and the observational curve (the 
solid line) are almost indistinguishable. 



This work was supported in part by Grant-in- Aid of Scientific Research of the Ministry of Education, 
Culture, and Sports, No.ll640274 and 09NP0801. 



''This might be possible if the statistical error approaches the systematic error in, for example, the SNAP 
project. 
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Fig. 1. — 54 Typela supernovae data from Supernova Cosmology Project (Perlmutter et al. 1999) in 
(i?, /) plane. / — + z and i? is a new independent variable defined by i? = {Ho/c){dL{z)/{l + z)), 

where Hq and dL{z) are the present Hubble parameter and the luminosity distance, respectively. The 

dotted line shows the flat dust case f — 1 — 0.5i?. The solid line is given by the likelihood analysis; 
/ = 1/V^, u = 1 + R + 0.084i?2 + 0.343i?3 + 0.360i?'' with = 47.72 for 54-4 d.o.f. The long and short 
dashed lines correspond to ±lcr values of U2,U3 and U4- 




0.5 1 

R 

Fig. 2. — 38 Simulated Typela supernovae data in {R,f) plane. See text how this data is 
made. The dotted line shows the flat dust case / = 1 — 0.5-R. The solid line is the result of 
likelihood analysis of / = 1/Vl + R + 0.2240i?2 + 0.1505i?3 + 0.0585i?4 + 0.0187i?5 + 0.0095i?6 + 0.0053i?7 
with = 36.16 for 38-7 d.o.f.. The dashed line is the theoretical curve of / = 

1/Vl + R + 0.225i?2 + 0.15i?3 + 0.054375i?4 + 0.016875i?5 + 0.0068906i?6 + 0.002257?^. Note that error 
bars are extended by a factor 100. 



